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Abstract 

The ROCket Combustor interactive Design 
(ROCCID) Methodology is a newly developed, interac- 
tive computer code for the design and analysis of a 
liquid propellant rocket combustion chamber. The 
application of ROCCID to design a liquid rocket com- 
bustion chamber is illustrated. Designs for a 
50 000 lbf thrust and 1250 psi chamber pressure com- 
bustor using Liquid Oxygen (L0X)/RP-1 propellants 
are developed and evaluated. Tradeoffs between key 
design parameters affecting combustor performance 
and stability are examined. Predicted performance 
and combustion stability margin for these designs 
are provided as a function of the combustor operat- 
ing mixture ratio and chamber pressure. 

Introduction 

An analytical design methodology to perform 
design tradeoffs leading to a high performing and 
stable combustion device for liquid propellant 
rocket engines has been recently developed. This 
design methodology, ROC ket Combustor interactive 
Design (ROCCID), contains previously developed and 
available analysis models for characterizing the 
effects of critical design and operating parameters 
on the performance and combustion/stability of liq- 
uid propellant combustors. Recently, ROCCID has 
been released to a limited number of users in a 
"Beta Test" (selected users debug and evaluate soft- 
ware before general release). 

In addition to the Beta Test evaluation, a sep- 
arate effort to validate ROCCID is underway. The 
objective of this validation effort is to use ROCCID 
to identify critical design and operating parameters 
that have a first order influence on combustion sta- 
bility and performance, to evaluate combustor 
designs created by ROCCID, and to demonstrate 
through hot-fire testing the extent to which ROCCID 
can be used as an "a priori" predictive tool. Some 
general rules were established to evaluate the hard- 
ware design for validation. First, the selected 
design must have the capability of providing for 
both stable and unstable operation through changes 
in operating conditions (chamber pressure and/or 
mixture ratio) as predicted by ROCCID. Ideally, 
this change between stable and unstable operation 
should be achievable with and without acoustic damp- 
ing devices. Second, the designs must be suffi- 
ciently different from designs tested previously so 
that a clear validation of ROCCID as an "a priori 
predictive tool is achieved. 
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ROCCID Methodology Overview 
ROCCID contains three main components: 

(1) An Interactive Front End (IFE) that pro- 
vides guidance to the user for input setup, input 
and output control, and the generation and mainte- 
nance of library files for replay, propellant prop- 
erties, and combustion gas properties. 

(2) A point analysis option that provides per- 
formance and combustion stability analysis of exist- 
ing combustor designs. 

(3) A point design option that creates the 
essential combustor design features for a high per- 
formance and stable rocket engine from specified 
design requirements. 

Both the point analysis and point design 
options use a library of performance and combustion 
analysis models that were selected from existing 
appropriate codes. These analysis models are con- 
tained within ROCCID in a modular format . This 
permits the user to access specific models for a 
specialized sub-analysis or to use two or more 
models that perform similar functions to define and 
resolve uncertainties in that particular area of the 
analysis. Modular construction also permits easier 
and less costly upgrades to the methodology as new 
analysis models are developed. 

The IFE is a menu driven preprocessor con- 
structed using an extensive library of interactive 
subroutines. Each input character is checked in 
the IFE for validity. Warnings are displayed when 
input errors are sensed. Replay files that contain 
all case input are. created and maintained. These 
files can be edited and used as input for a subse- 
quent session. The user may repeatedly alter the 
input until the desired result is achieved. 

Required propellant properties for both the injected 
fuel and oxidizer and their resulting combustion gas 
are internally generated and maintained in files. 

Upon completion of the input deck, the user 
may execute either the point analysis or point 
design options. These options contain many interac- 
tive decision points for the user. Upon completion 
of an analysis or user termination, the run stream 
returns to the IFE to monitor the results, record, 
print pertinent information, and prepare input for 
the next run. 

The performance of the combustor is defined by 
the energy release efficiency. This accounts for 
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combustion efficiency limitations resulting from 
incomplete propellant vaporization and/or mixing. 

The energy release efficiency is calculated using 
the JANNAF simplified performance procedure 3 and 
the propellant vaporization and mixing limitations 
from the steady state combustion analysis. An 
input file for use in the TDK/BLM performance analy- 
sis of the complete engine is also generated. 

A steady state combustion analysis including 
propellant atomization, vaporization, and mixing 
supplies key input to the performance and stability 
analyses. Four models/correlations for propellant 
drop size are included for standard injector ele- 
ments including showerhead, doublet and triplet 
impinging elements, and shear and swirl coaxial ele- 
ments, Drop sizes from all applicable correlations 
are calculated and displayed for comparison. A 
user may select any of the calculated values for 
performance and stability analyses or provide 
another estimated value. Propellant vaporization 
is calculated using the generalized length correla- 
tion. 2 Propellant mixing is based on the use of a 
unielement mixing efficiency value determined previ- 
ously from cold flow measurements and adjusted for 
interelement mixing effects. 

The combustion stability analysis can be made 
with a large array of models used to calculate the 
chamber admittance and the burning and injector 
response magnitudes. These models provide the capa- 
bility to estimate combustion stability margin for 
all common types of combustor instabilities includ- 
ing chugging and chamber acoustic coupled (high fre- 
quency) modes. The effects of damping devices such 
as acoustic cavities and baffles are also considered 
by these models. 

The design requirements for combustor cooling 
are established using techniques outside of ROCCID. 
These requirements may include estimates of fuel 
film cooling required for chamber and baffle walls, 
dump cooling of baffle tips, and bulk temperature 
increases resulting from regenerative cooling of 
the nozzle, chamber, and resonator/baffle compo- 
nents. This information is used to calculate the 
propellant injection temperatures, in jec t ion or i f ice 
requirements, and the local flow injection mixture 
ratios . 

ROCCID displays the stability analysis results 
in terms of the calculated growth coefficient (X) 
for the particular acoustic mode of concern. This 
growth coefficient represents the amount of amplifi- 
cation required by the chamber to achieve the condi- 
tion where the driving required to support the 
waveform exactly equals the driving response present 
within the system. Thus a value of X - 0 repre- 
sents a neutral stability condition, X > 0 repre- 
sents an unstable condition, and X < 0 represents 
a stable operating condition. 

The growth coefficient, in terms of the tradi- 
tional combustion instability transfer function for- 
mat, is illustrated in Fig. 1. Three typical 
situations are illustrated in this example. 

Figure 1(a) illustrates results of a tradi t ional 
stabi li ty analysis where the shaded region formed 
by overlapped chamber response and combustion 
response curves indicate the potential for an insta- 
bility. However, this traditional analysis does not 
include a survey of chamber admittances that would 
lead to a growth coefficient value which would pro- 
vide zero stability margin, i.e. r chamber response 


and combustion response curves being tangent at some 
point. The growth coefficient at this zero margin 
point provides a quantitative assessment of stabil- 
ity. Figures 1(b) and (c) illustrate unstable and 
stable quantitative predictions, respectively, 
wherethe growth coefficient is used to define the 
relative stability based on the rate of pressure 
applitude growth (X > 0) or decay (X < 0). The fol- 
lowing is the exponential relationship used to 
define X: 

A/A 0 = 

where A/A 0 is the pressure amplitude ratio change 
over the period, At (sec). 

The resultant growth coefficient can be related 
to the damp time required for dynamic stability as 
defined in CP1A 247, by defining a growth coeffi- 
cient (X£p{£ 247) required to achieve a 10:1 over- 
pressure damping within the damp time specified from 
CPIA 247. That is: 

’CPIA 247 - I" (0 1)/tdajnp, CPIA 247 

where 

tdamp, CPIA 247 - l .250//f (sec), 

f - frequency of the unstable oscillation 

reducing 

*CPIA 247 » -1*842 /f (1/sec) 

Mote that ^CP[\ 247 w ill always be negative 
indicating a damped system and will be a function 
of the frequency (mode) of the damped oscillation. 
Therefore, if the growth coefficient is equal to or 
less than X(]p[A 247 the system is dynamically sta- 
ble. Conversely, if the growth coefficient is 
greater than X^pj^ 247 the system is not dynami- 
cally stable. 

Validation Design Requirements 

The creation of a combustor design consists of 
an evolutionary process where design requirements 
and operational goals are used to establish design 
parameters for the injector and combustion chamber. 
Iterations to the set of critical design parameters 
are performed to meet performance goals while pro- 
viding for stable combustion. ROCCID has been used 
to examine the sensitivity and influence of several 
injector and combustor design parameters on the pre- 
dicted combustor performance and combustion 
stabi 1 i ty . 

The propellant combination of L0X/RP-1 has been 
selected for study using ROCCID and for eventual 
hot-fire testing to validate ROCCID capabilities. 

This propellant was selected because historically, 
it has proven to be a difficult propellant to pro- 
vide both stable combustion and high performance. 

Nominal design propellant injection tempera- 
tures of 174.4 °R for liquid oxygen and ambient con- 
dition (515 °R) for RP-1 have been selected. These 
conditions represent normal propellant storage tem- 
peratures within the Aerojet test facility and 
therefore are the most cost effective test 
condi t ions . 
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A nominal design mixture ratio of 2.8 has been 
selected since this approximates an optimum perform- 
ance design from a thrust chamber performance view- 
point. Mote that this mixture ratio value has been 
selected for design purposes only. The effects of 
mixture ratio variation on a fixed design have been 
evaluated using ROCCID and will be a key operating 
parameter for validation testing. 

Thrust chamber size is a key design parameter, 
particularly from a combustion stability standpoint. 
The chamber diameter has a direct effect on the res- 
onant frequency of acoustic modes within the cham- 
ber. Large chambers appropriate for booster engine 
application (Dc - 17.5 - 44 in.), have lower reso- 
nant frequencies that are more likely to couple with 
typical combustion responses. On the other hand, 
subscale chamber diameters of 5.5 and 7.68 in. have 
been successfully used in the past for combustion 
stability investigations . 4 For this study, a nom- 
inal baseline chamber diameter of 7.68 in. has been 
selected since it is a proven '’subscale ’ design and 
residual hardware is available from the test program 
in Ref. 5 , 

A contraction ratio of 2:1 that corresponds to 
a throat diameter of 5.43 in. and a nominal design 
chamber pressure of 1250 psi were selected to pro- 
vide a large operating chamber pressure variation 
within the existing hardware and test fac i 1 i ty capa- 
bilities. The operational envelope of the valida- 
tion hardware for the Aerojet E-4 test facility is 
shown in Fig. 2. A variation in mixture ratio from 
approximately 1 to 10 is possible at the nominal 
chamber pressure of 1250 psi. A maximum chamber 
pressure of approximately 1800 psi is possible at 
the nominal mixture ratio. The mixture ratio 
excursion range diminishes as chamber pressure is 
increased due to the facility run tank pressure 
1 imi ts . 

Impinging injector designs were selected for 
this evaluation since they are most appropriate for 
liquid-liquid injection. Specifically, the like 
doublet and 0F0 triplet elements, both of which are 
included in ROCCID, were evaluated. The 0F0 triplet 
element was selected over the F0F since the propel- 
lant density and mixture ratio will result in equal 
orifice diameters with an 0F0 element but not with 
the F0F triplet. From experience, equal orifice 
diameters are desirable with a triplet element to 
provide optimum mixing and atomization. 

A summary of the selected design parameters 
for this study is presented in Table 1. Design 
parameters from two recent combustion stability 
investigations in Refs. 4 and 5, are also provided 
for comparison purposes. 

Downselect Criteria 

Criteria were established to select which 
injector to use to validate ROCCID. Because of lim- 
ited resources, one and at best two different thrust 
chamber designs could be selected for design and 
testing. The criteria were established on the basis 
of providing the best test of the capabilities of 
ROCCID. The criteria are: 

(1) Largest negative and positive growth coef- 
ficients for the most stable and unstable modes 
without acoustic damping. 


(2) Highest sensitivity of combustion stabil- 
ity To chamber pressure and mixture ratio. 

(3) Greatest design change from existing data 

base . 

(4) Capable of stabilization at nominal operat- 
ing conditions using an acoustic cavity. 

(5) Highest confidence in modeling drop size 
and mixing efficiency. 

(6) Ease of incorporating design features into 
exist ing hardware . 

Selection of the Injector 

Three injector designs, Fine triplet, coarse 
triplet, and like-on-like doublet, were selected to 
be evaluated by ROCCID for detail design and test- 
ing. Initially, the point design option of ROCCID 
was used to size the injectors for the nominal oper- 
ating conditions. Thereafter, the analysis portion 
of ROCCID was used to perform trade studies on the 
initial designs. The experience gained from the 
experimental results of Ref. 4 was also applied 
to the analysis. Also from past experience, certain 
model combinations in ROCCID were used because they 
produced good results. H I F I was used to analyze 
the chamber response. Smi th-Reardon and N/x cor- 
relations were used to analyze the burning response. 
INJ was used to analyze the injector response. 

The injector orifice sizes evaluated were 
0.090 in. for the fine triplet, 0.159 in. for the 
coarse triplet, and 0.100 in. L0X and 0.065 in. 
fuel for the like-on-like doublet. The predicted 
performance efficiencies For all three injectors 
were greater than 97 percent so performance did not 
become a critical concern in the selection process. 

An extensive number of stability predictions 
were generated and only a portion will be discussed 
in the following section. References 6 and 7 should 
be examined for more details. 

The chug stability of the three injectors was 
analyzed. The fine triplet and the doublet have 
low chug thresholds at about 300 psi. The coarse 
triplet has a more undesirable chug threshold at 
about 600 psi . Considering that the nominal operat- 
ing pressure is 1250 psi, this higher chug pressure 
limit for the coarse triplet does not allow for much 
leeway in throttling the engine to find stable and 
unstable high frequency stability regions. 

The high frequency stability characteristics 
of the three injectors were also analyzed. The 
coarse triplet was found to be the most stable 
injector design (X less than -37 in all modes) 
without acoustic damping devices. From a flight 
hardware designer's point of view, this would be an 
excellent design. However, since our validation 
cri teria^ requi re both ""stable and unstable operation 
for verification, the coarse triplet injector design 
is not sufficient for ROCCID validation. If 
resources are available, this design may turn out 
to be a good choice for a secondary, very stable, 
alternate test series to demonstrate ROCCID' s capa- 
bility to stabilize the thrust chamber through 
injector design changes (i.e., combustion response). 
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The doublet was found to be stable except for 
the first tangential mode. The doublet looks very 
promising because it has a region of operation 
where it can be driven stable and unstable by chang- 
ing the mixture ratio as shown in Fig. 3. The dou- 
blet can be stabilized using acoustic damping 
devices which will allow validation of the models 
for damping devices. The major drawback of the dou- 
blet is that it is very similar to the design used 
in Ref. 5. Since this is a validation of the pre- 
dictive capability of the models in ROCCID, using 
an injector so close to the anchored hardware makes 
this injector a secondary choice. 

The fine triplet design was found to be unsta- 
ble in several different modes with varying 
operating conditions. Figures 4 to 7 show ROCCID 
predictions of stability changes over the operating 
map. The first longitudinal (1L) mode is stable 
until lower chamber pressures are reached. The 
first tangential (IT), second tangential (2T), and 
first radial (1R) modes are unstable over wide 
changes in operating conditions. The orifice dia- 
meter was changed from 0.09 to 0.10 in. to see what 
the effect of orifice diameter had on stability. 

It did not improve stability margin sufficiently to 
provide any significant advantages. Damping devices 
were found to be very effective in damping out the 
unstable modes. Also, this design is very different 
from the hardware described in Reference 5 in terras 
of number of elements and orifice sizes. 

Table 2 provides a rating of the candidate 
injector designs. It is obvious that the fine tri- 
plet injector is the best choice to validate the 
models contained within ROCCID. The fine triplet 
injector stability is sensitive to chamber pressure 
and mixture ratio variations which will allow for 
testing the sensitivity and predictive capabilities 
of the models. The fine triplet injector is suffi- 
ciently different from previously tested injectors 
to allow a true "a priori' 1 prediction. As a result, 
this injector was selected as the test hardware. 

Growth Coefficient Uncertainty 

The sensitivity of the stability predictions 
(in terms of the growth coefficient) to the pres- 
sure interaction index, N, and sensitive time lag, 
x, were evaluated over a range of N and x val- 
ues. The ranges of N and x values were selected 
based on estimates from three sources. 

(1) The empirical variations in these parame- 
ters shown in Refs. 8 and 9, 

(2) The calculated variation resulting from 
different droplet combustion mixture ratios and gas 
temperatures using the CRP model. ^ 

(3) Variations in x based on an estimated 1 
percent uncertainty in energy release efficiency and 
relating that uncertainty to a change in vaporiza- 
tion time lag. 

The three techniques yielded similar results. 
Based on the Ref. 8 data the uncertainty in N was 
estimated to be +26 percent (0.80+0.20) and the 
uncertainty for x was +43 percent, -31 percent 
(0.7xl0- 4 sec to 1.79xl0- 4 sec) 8 (Figs. 6(c) 
and (d)). 

The CRP model was run for assumed fuel droplet 
surface mixture ratios of 0.7 to 2.4 and gas temper 


atures of 2500 to 6500 °R at the droplet surface. 
These analyses resulted in N values of 0.45 to 
1.0 and x values of 0.94xl0- 4 sec to 2.5xi0- 4 sec. 
similar to the ranges shown from Ref. 8 but slightly 
larger . 

Assuming a 1 percent performance uncertainty 
(98+1 percent) caused by a fuel droplet size uncer- 
tainty would yield a x uncertainty of approxi- 
mately 20 percent, again in the general range 
observed from the Smi th-Reardon data. 

The ROCCID model was run with the above varia- 
tions in N and x f rom/Srai th-Reardon at the nominal 
chamber pressure and mixture ratio (P c - 1250 psi, 
0/F - 2.8). The results, as shown in Fig. 8 indi- 
cate large variations in predicted growth coeffi- 
cient, up to approximately +600 sec -1 for the 2T 
mode if both N and x variations are considered. 
However, implicit within ROCCID is a variation in x 
due to engine operating conditions. Therefore, the 
data scatter leading to uncertainties in x from 
Refs. 8 and 9 may be due to differences in the oper- 
ating condition of the engines being evaluated. If 
only the uncertainty in N is considered, the 
growth coefficient variations are reduced to 
+25 sec' 1 for the 1L mode, +350 sec' 1 for the IT 
mode, and +500 sec' 1 for the 2T mode. Based on the 
statistical scatter in the existing empirical corre- 
lations of N and the present state of the models, 
these large bands of uncertainty associated with 
the predictions are expected. 

It should be noted that there is not an equal 
expectation of N values within the published un- 
certainty band of +26 percent. The nominal value 
is the most likely with the extremes, +26 percent, 
probably representing a +3a standard deviation. If 
this is the case, then a statistical analysis of 
experimental data going into ROCCID could be per- 
formed to estimate the confidence level of the ROC- 
CID predictions for calculated values of growth 
coefficient. A more detailed analytical and statis- 
tical evaluation of the data is recommended to allow 
a better understanding of the reasons for and conse- 
quences of the uncertainties in H and x . 

Validation Test ing 

Test planning is proceeding at this time. 

Three test points without damping devices have been 
determined. Looking at Figs. 4 to 7, one can see 
that at chamber pressure/mixture ratio combinations 
of 625 psi/7.50, 1250 psi/1.15, and 1250 psi/7.50 
stability should be achieved in all modes. These 
points will be tested to determine the model's capa- 
bility of predicting undamped stable points. Sev- 
eral of the other points shown in Figs. 4 to 7 will 
be tested both with and without damping devices to 
determine the effectiveness of the damping device 
and predictive capability of the models. 

Conclusions 

It is evident from the extent of trade studies per- 
formed in a very brief period that ROCCID is a pow- 
erful tool. The ease with which an engine designer 
can make and analyze changes can potentially save 
many manhours. ROCCID is valuable because it docu- 
ments how the designer arrived at a particular 
design. In the future, validation of the models in 
ROCCID can lead to upgrading insufficient models and 
determining the better analysis models. From this 
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analysis, it is evident that the band of uncertainty 
associated with the models that industry relies on 
to design flight hardware is too broad and requires 
further study. The planned test program based on 
the trade studies will determine the "a priori'' pre- 
dictive capabilities of the models. 
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TABLE 1 . - COMPARISON OF SELECTED NOMINAL DESIGN PARAMETERS 


Combustion 

stability investigation 



NAS 

3-25556 a 

F046U-85-C-0100 b 

NAS 3-24612 c 

Propellant combination 


L0X/RP-1 

L0X/RP-1 

LOX/CH 4 





„ - 1 .6 

Mixture ratio 


2 , 8 iTBD 

2 .9±0 . 3 

3,0 +0.2 

Injection temperature 
L0X °R 


174 

174 

174 

Fuel °R 


515 

530±20 

530 to 438 

Chamber pressure, psia 


1250±TBD 

2000 -S 

+100 
^ uuu _500 

Chamber diameter, in. 


7.68 

7.68 

5.66 

Throat diameter, in. 


5.43 

4.43 

3.31 

Contraction ratio 


2:1 

3:1 

2.92:1 

Chamber length, in. 


TBD 

8.5,13,20 

14 

Nominal thrust level, 
lbf 


48 000 

57 000 

30 000 

Total flowrate' 
lbm/sec 


150±TBD 

180 -90 

190 -| 5 7 

Injection element type 

0F0 

triplet 

0F0 triplet 

Shear coaxial 

1 ike 

doub 1 e t 

1 ike doublet 



a L0X/Hydrocarbon Rocket Engine Analytical Design Methodology Development and 
Validation 

^Oxygen /Hydrocarbon Injector Characterization^ 
c L0X/Hydrocarbon Combustion Instability Investigation 4 
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TABLE 2. - INJECTORS COMPARED ACAINST RATING CRITERIA 


Rating Cr i ter ia 

Injector configuration 

Fine 0F0 
triplet 

Coarse 0F0 
triplet 

Li ke-on-1 i ke ; 
one doublet 

Largest negative and positive 
growth coefficients for least 
damped mode without acoustic 
damping 

’ ’ .'tips t unstable 

v ’ ’ tiqsf ; s 

Unstable 

Best sensitivity of combustion 
stability to P c and 0/F 

Significant variation 
; ; . 2T and ’ it ; s t able 

aeration with 0/F and 

1 p c :?■ • 

Much wi lder 
variation in P c 
and 0/F required 
to drive unstable 

: 5ppd sensitivity 
^0/P for IT 

Greatest change from existing 
data base 

Ne , 39 to XOS 
D 0 , 0.125 to 0.090 

Ne, 39 to 48 
D 0 . 6 .125 to 0tS9 

Ne, 105 to 123 
0 o , 0.059 to 0.065 

Capable of stabilization at 
nominal operating conditions 
using an acoustic cavity 

Three types of 
cavities evaluated ' 

Stable wi thout an 
acoust ic cavi ty 

fT mono t ime /’ 

: cavity judged to 
: be .sufficient 

Best confidence in modeling 
drop size and mixing 
ef f iciency 

• r® reasonable ; 
C good 

r™ reasonable; 

En. good 1 mi . 

r m reasonable; 
E m reasonable 

Ease of incorporating design 
features into existing hard- 
ware 

New faceplate only 

New faceplate and 
injector core 

New faceplate and 
injector core 

Preferred rating 

6 of 6 

3 of 6 

2 of 6 


Indicates prefered characteristic 
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Fig. 3— Like doublet injector stability prediction, chamber pressure = 1250 psia; 
length - 17 in.; no damping. 
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Fig. 4, — Fine triplet 1 L mode stability map showing values of X at various points as 
well as point predicted to be extremely stable (S). 



Growth 

coefficient, 

X 

= 0 

Test facility limits 
Nominal chamber pressure 
Nominal mixture ratio 


Chamber pressure, P c , psia 


Fig. 5.— Fine triplet IT mode stability map showing values of X at various points as 
welt as points predicted to be extremely stable (S). 
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Fig.6 — Fine triplet 2T mode stability map showing values of X at various points as 
well as points predicted to be extremely stable (S). 
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Fig. 7 —Fine triplet 1 R mode stability map showing values of X at various points as 
well as points predicted to be extremely stable (S). 
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Fig. 8. — Variation of growth coefficient X with N and x. 
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